What is claimed is: 
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1. A method of scaling up permeabilities associated with a fine-scale grid of cells 
representative of a porous medium to permeabilities associated with an 

5 unstructured coarse-scale grid of cells representative of the porous medium, 

comprising the steps of: 

(a) generating an areally unstructured, Voronoi, computational grid using 
the coarse-scale grid as the genesis of the computational grid, the cells of 
the computational grid being smaller than the cells of the coarse-scale 

10 grid and each cell of said computational grid and said coarse-scale grid 

having a node; 

(b) populating the computational grid with permeabilities associated with 
the fine-scale grid; 

(c) developing flow equations for the computational grid, solving the flow 
1 5 equations, and computing inter-node fluxes and pressure gradients for 

the computational grid; 

(d) using the fluxes and pressure gradients computed in step (c) to calculate 
inter-node average fluxes and average pressure gradients associated with 
the coarse-scale grid; and 

20 (e) calculating upscaled permeabilities associated with the coarse-scale grid 

using the average fluxes and average pressure gradients calculated in 
step (d). 

2. The method of claim 1 wherein inter-node connections of the computational 
25 grid cells are parallel to inter-node connections of the coarse-scale grid. 

3. The method of claim 1 wherein cells of the computational grid are 
approximately the same size as the cells of the fine-scale grid. 

30 4. The method of claim 1 wherein cells of the computational grid are smaller 
than the cells of the fine-scale grid. 
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5. The method of claim 1 wherein the permeabilities populated in step (b) are 
assigned to nodes of the computational grid, the permeability assigned to a 
given node of the computational grid corresponding to a predetermined 
5 permeability of a cell of the fine-scale grid that would contain the given node's 

location if the computational grid was superimposed on the fine-scale grid. 



6. The method of claim 5 further comprises calculating the permeability for a 
given inter-node connection of the computational grid by harmonically 
10 averaging the permeabilities at the two nodes forming the given inter-node 

connection. 



7. The method of claim 1 wherein the permeabilities populated in step (b) are 
assigned to midpoints of a given inter-node connection of the computational 
1 5 grid, the permeability assigned to the given inter-node connection 

corresponding to a predetermined permeability of a cell of the fine-scale grid 
that would contain the mid-point of such inter-node connection if the 
computational grid was superimposed on the fine-scale grid. 

20 8. The method of claim 1 wherein the calculation of the inter-node average 

fluxes and average pressure gradients in step (d) is determined by using only 
the fluxes and pressure gradients computed for the inter-node connections of 
the computational grid that fall within a predetermined sub-domain of the 
computational node. 

25 

9. The method of claim 1 wherein the calculation of the average fluxes and 

average pressure gradients in step (d) for a given inter-node connection of the 
coarse-scale grid is determined using only fluxes and pressure gradients 
computed in step (c) for the inter-node connections of the computational grid 
30 that are parallel to the given inter-node connection. 
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10. The method of claim 1 wherein the permeabilities calculated in step (e) are 
determined for a given node of the coarse-scale grid. 

1 1 . The method of claim 1 wherein the permeabilities calculated in step (e) are 
5 determined for a given inter-node connection of the coarse-scale grid by 

calculating the ratio of the average flux to the average pressure gradient 
computed in step (d) for the given inter-node connection. 

12. The method of claim 1 wherein cells of the fine-scale grid are structured. 

10 

13. The method of claim 1 wherein the coarse-scale grid is a PEBI grid. 

14. The method of claim 1 wherein both the coarse-scale grid and the 
computational grid are PEBI grids. 

15 

15. The method of claim 1 wherein inter-node connections of the coarse grid forms 
Delaunay triangles and the computational grid generated in step (a) contains 
similar, smaller Delaunay triangles equal in number to n 2 times the number of 
Delaunay triangles of the large-scale grid, where n is a predetermined integer 

20 refinement ratio used to generate the computational grid. 

16. The method of claim 1 wherein all cells are three-dimensional. 

17. The method of claim 16 wherein the coarse-scale grid and the computational 
25 grid are both unstructured areally and structured vertically. ^ 

18. The method of claim 1 further comprises determining inter-node connection 
transmissibilities of the coarse-scale grid using permeabilities calculated in 
step (e). 

30 
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19. The method of claim 1 wherein the flow equations of step (c) are single-phase 
and steady-state. 



20. A method for estimating permeability of each cell of a first grid having a 

5 multiplicity of cells of a subterranean geologic domain using a predetermined 

permeability for each cell of a second grid representative of the domain, said 
second grid containing a larger number of cells than the first grid, the method 
comprising: 

(a) constructing an unstructured, third grid representative of the domain 
10 comprising approximately the same or greater number of cells than the 

second grid, each cell of the first, second, and third grids having a node 
and each link between two nodes of adjacent cells being an inter-node 
connection, substantially all of the inter-node connections of the third 
grid being parallel to the inter-node connections of the first grid; 
1 5 (b) for each node of the third grid, assigning a permeability corresponding to 

the permeability of a cell of the second grid that contains the node 
location of the third grid; 
(c) developing a single-phase, steady-state pressure equation for each cell of 
the third grid system; 
20 (d) solving the pressure equations and computing fluxes and pressure 

gradients for all inter-node connections of the third grid; 

(e) computing an estimated permeability for a given connection of the first 
grid using inter-node connections of the third grid; and 

(f) repeating step (e) for all connections of the first grid. 

25 

21 . The method of claim 20 further comprises computing the permeability in 
step (e) by the additional steps of determining average fluxes and average 
pressure gradients over sub-domains associated with a given grid inter-node 
connection of the first grid and calculating a ratio of the average flux to the 

30 average pressure gradient, thereby obtaining the permeability for the given 

inter-node connection of the first grid. 
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22. The method of claim 20 wherein the permeability computation of step (e) uses 
only inter-node connections of the third grid that are parallel to inter-node 
connections of the first grid. 

23. The method of claim 20 wherein all grid cells are three-dimensional. 

24. The method of claim 23 wherein the second grid and the third grid are each 
unstructured areally and structured vertically. 

25. A method for estimating permeabilities associated with cells of a large-scale 
grid representative of fluid flow in a porous medium using predetermined 
permeabilities associated with cells of a small-scale grid also representative of 
fluid flow in the porous medium, each cell of the large-scale grid having a 
node and each node being linked to adjacent nodes to form inter-node 
connections and such connections forming Delaunay triangles, comprising the 
steps of: 

(a) constructing a computational grid by dividing each Delaunay triangle of 
the large-scale grid into a multiplicity of similar, smaller Delaunay 
triangles, the sides of such smaller Delaunay triangles being inter-node 
connections of the computational grid and the inter-node connections of 
the large-scale grid and the computational grid being aligned with each 
other; 

(b) assigning permeabilities to the computational grid corresponding to the 
predetermined permeabilities of the small-scale grid; 

(c) developing a single-phase, steady-state pressure equation for each cell of 
the computational grid, solving the pressure equations, and computing 
fluxes and pressure gradients for all inter-node connections of the 
computational grid; 
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(d) using the fluxes and pressure gradients computed in step (c) to calculate 
an average flux and an average pressure gradient for each inter-node 
connection of the large-scale grid; and 

(e) calculating a permeability associated with a given inter-node connection 
of the large-scale grid using the average flux and average pressure 
gradient calculated in step (d). 

26. The method of claim 25 wherein the number of nodes of the computational 
grid are approximately the same as the number of cells of the small-scale grid. 

27. The method of claim 25 wherein the number of nodes of the computational 
grid are smaller than the number of cells of the small-scale grid. 

28. The method of claim 25 wherein the number of similar, smaller Delaunay 
triangles is n 2 times the number of Delaunay triangles of the large-scale grid, 
where n is a predetermined integer refinement ratio. 
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